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Abstract 

The gauge independence of the dynamical fermion mass generated through chiral 
symmetry breaking in QED in a strong, constant external magnetic field is critically 
examined. We show that the bare vertex approximation, in which the vertex cor- 
rections are ignored, is a consistent truncation of the Schwinger-Dyson equations in 
the lowest Landau level approximation. The dynamical fermion mass, obtained as 
the solution of the truncated Schwinger-Dyson equations evaluated on the fermion 
mass shell, is shown to be manifestly gauge independent. By establishing a di- 
rect correspondence between the truncated Schwinger-Dyson equations and the 2PI 
(two-particle-irreducible) effective action truncated at the lowest nontrivial order in 
the loop expansion as well as in the 1/iVy expansion (Nf is the number of fermion 
flavors), we argue that in a strong magnetic field the dynamical fermion mass can 
be reliably calculated in the bare vertex approximation. 
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1 Introduction 



The study of chiral symmetry breaking has a long history, beginning with 
the seminal work of Nambu and Jona-Lasinio [1-3]. Because of its simpler 
structure, quantum electrodynamics (QED) has widely been used to study 
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chiral symmetry breaking in gauge theories [4-12]. Being inherently a nonper- 
turbative phenomenon, the dynamical generation of a fermion mass in gauge 
theories is usually studied with the help of the Schwinger-Dyson equations 
truncated in certain schemes. An often used truncation scheme is the rainbow 
(or quenched ladder) approximation (for a review see Ref. [13]). It is gener- 
ally difficult to find a consistent truncation scheme such that the resulting 
dynamical fermion mass is gauge (fixing) independent. A common practice 
has been to forgo gauge independence in favor of obtaining a nonperturba- 
tive solution. Nevertheless, the demonstration of the gauge independence of 
physical quantities is of fundamental importance in gauge theories [14,15]. In 
perturbation theory the usual connection between the order of the loop ex- 
pansion and powers of the gauge coupling guarantees that physical quantities 
calculated to a given order in the loop expansion are gauge independent. In 
other words, as is well known, loop expansion is a consistent expansion scheme 
in perturbation theory. Such a connection however is lost in nonperturbative 
studies in which contributions of a given order in the gauge coupling arise 
from every order in the loop expansion. As a result, establishing the gauge 
independence of physical quantities obtained in a nonperturbative calculation 
(e.g., the Schwinger-Dyson equations) is a highly nontrivial problem. This is 
because an infinite subset of diagrams arising from every order in the loop 
expansion has to be resummed consistently, thus leading to potential gauge 
dependence of physical quantities whenever not all relevant diagrams are ac- 
counted for. Therefore, we emphasize that in gauge theories no truncation 
schemes of the Schwinger-Dyson equations should be considered consistent 
unless the gauge independence of physical quantities calculated therein is un- 
equivocally demonstrated. 

The problem of chiral symmetry breaking in QED in a strong, constant exter- 
nal magnetic field has received a lot of attention over the past decade [16-29] . 
It was originally motivated by a proposal [30-33] to explain the correlated 
e + e~ peaks observed in heavy ion collision experiments [34,35]. The proposal 
posits that the e + e~ peaks are produced from the decay of a bound e + e~ 
system formed in a new metastable phase of QED with broken chiral sym- 
metry, which is induced by the strong electromagnetic field present in the 
neighborhood of the colliding heavy ions. Even though the e + e~ events were 
not observed in a later experiment [36], it is still an interesting question to ask 
whether background fields can induce chiral symmetry breaking. The result 
of the investigation suggests that it may be relevant for studying the elec- 
troweak phase transition in the early universe when a strong magnetic field 
was present [18,19]. The methodology developed also found interesting ap- 
plications in the study of color superconductivity in high density QCD [37], 
superconductivity and magnetic field induced phase transitions in condensed 
matter systems [38-42], as well as the properties of strongly magnetized as- 
trophysical media [43,44]. 
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In a recent article [45], we have found a consistent truncation scheme for 
the case of chiral symmetry breaking in a strong, constant external magnetic 
field. We critically examined the consistency and gauge independence of the 
bare vertex approximation that has been extensively used in truncating the 
Schwinger-Dyson equations to calculate the dynamical fermion mass gener- 
ated through chiral symmetry breaking in weakly coupled QED in a strong 
magnetic field. In particular, we showed that the bare vertex approximation, 
in which the vertex corrections are ignored, is a consistent truncation of the 
Schwinger-Dyson equations in the lowest Landau level approximation. The 
dynamical fermion mass, obtained as the solution of the truncated Schwinger- 
Dyson equations evaluated on the fermion mass shell, is manifestly gauge 
independent. This novel gauge independent approach poses a serious question 
on the validity of the results and conclusions obtained in earlier studies [18- 
20,22,24-29]. Most importantly, it leads to a consistent, nonperturbative cal- 
culation of the dynamical fermion mass to leading order in the gauge coupling, 
and allows one to identify the infinite subset of diagrams that contribute to 
chiral symmetry breaking in a strong magnetic field. 

The purpose of this article is twofold. On the one hand, we present the tech- 
nical details of the results outlined in Ref. [45]. To begin with, we critically 
review the properties of the so-called E p functions, first introduced by Ri- 
tus [46,47] in the studies of QED in constant external electromagnetic fields. 
In doing so, we correct a few mistakes found in the literature [20,22,23,46-48]. 
We also derive the Ward-Takahashi identity in the bare vertex approximation 
both within and beyond the lowest Landau level approximation. This allows 
us to (i) correct a few oversights found in Ref. [23], (ii) verify the conclusion 
obtained therein on the Ward-Takahashi identity in the bare vertex approxi- 
mation and within the lowest Landau level approximation, and (iii) show that 
the Ward-Takahashi identity in the bare vertex approximation can be satisfied 
only within the lowest Landau level approximation. By utilizing the E p func- 
tion formalism and the Ward-Takahashi identity in the bare vertex approxima- 
tion, we prove that the bare vertex approximation is a consistent truncation of 
the Schwinger-Dyson equations in the lowest Landau level approximation. In 
particular, we verify that (i) the truncated vacuum polarization is transverse, 
(ii) the truncated fermion self-energy evaluated on the fermion mass shell is 
manifestly gauge independent. 

On the other hand, we provide a detailed comparison between our results and 
those obtained in the literature [18-20,22,24-29]. Based on the gauge inde- 
pendent analysis presented in this article, we argue that to a large extent the 
results and conclusions obtained in those earlier studies can be attributed to 
inconsistent truncation schemes, gauge dependent artifacts, or both. Along 
the way this detailed comparison also establishes several important aspects 
of the bare vertex approximation: (i) The assumption of a momentum inde- 
pendent fermion self-energy, a consequence of the Ward-Takahashi identity in 
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the bare vertex approximation, is reliable in the momentum region relevant to 
chiral symmetry breaking in a strong magnetic field, (ii) While the extension 
to a momentum dependent fermion self-energy violates the Ward-Takahashi 
identity in the bare vertex approximation, the resultant gauge dependence of 
the dynamical fermion mass is of higher order, (iii) In a strong magnetic field 
the dynamical fermion mass can be reliably calculated in the bare vertex ap- 
proximation and is gauge independent up to corrections of higher order, (iv) 
There exists a direct correspondence between the Schwinger-Dyson equations 
truncated in the bare vertex approximation and the nonperturbative approach 
based on the 2PI (two-particle-irreducible) effective action truncated at the 
lowest nontrivial order both in the loop expansion as well as in the 1/Nf 
expansion (Nf is the number of fermion flavors). 

The rest of this article is organized as follows. In Sec. 2 we critically review the 
properties of the Ritus E p functions. In Sec. 3 the Schwinger-Dyson equations 
truncated in the bare vertex approximation are discussed, and the Ward- 
Takahashi identity in the bare vertex approximation is derived both within and 
beyond the lowest Landau level approximation. In Sec. 4 we present the proof 
that the bare vertex approximation is a consistent truncation of the Schwinger- 
Dyson equations in the lowest Landau level approximation. The truncated on- 
shell Schwinger-Dyson equations are numerically solved for various numbers of 
fermion flavors, and an analytic fit for the dynamical fermion mass is obtained. 
In Sec. 5 we compare our results to those found in the literature. We also 
establish a direct contact with the 2PI effective action truncated at the lowest 
nontrivial order in the loop expansion as well as in the 1/Nf expansion. Finally, 
we present our conclusions in Sec. 6. In Appendix A we show that the vacuum 
current vanishes identically in a constant external magnetic field. Appendix B 
provides a summary of some exact relations derived from the 2PI effective 
action that are useful to the discussion in the main text. 



2 Formalism: Ritus E p functions 



In this section we critically review the properties of the E p functions, first 
introduced by Ritus in his seminal work [46,47]. They form a complete set 
of Dirac matrix- valued orthonormal functions and hence provide a convenient 
formalism in the studies of QED in the presence of a constant external mag- 
netic field [20,22,23,48]. In doing so, a few mistakes and confusions found in 
the literature [20,22,23,46-48] are corrected and clarified. 

We begin with the free field Dirac equation for a massless fermion in the 
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presence of a constant external magnetic field 1 

( 7 .7r)#r) = 0, (2.1) 

where 7r^ = — id^ — eA^ with being the external gauge field and ip is a 
Dirac spinor. Instead of solving Eq. (2.1) directly, we show that there exists a 
complete set of orthonormal functions, referred to in the literature as the Ritus 
E p functions, such that in the basis spanned by the E p functions Eq. (2.1) is 
rendered formally identical to the Dirac equation for a massless fermion in the 
absence of external fields. 

We will take the constant external magnetic field of strength H in the £3- 
direction with the corresponding vector potential given by 

^ = (0,0,^,0), (2.2) 

where \i = 0,1,2,3. It is straightforward to verify that the operators (7 • 
7r) 2 , E 3 = «7 1 7 2 and 7 5 constitute a maximal set of mutually commuting 
operators. 2 The E p functions are constructed in terms of the simultaneous 
eigenfunctions (eigenvectors) of these mutually commuting operators. 

The eigenfunction equation for (7 • 7r) 2 is given by 

-^■n) 2 <p p (x)=p 2 <p p (x), (2.3) 

where <f> p (x) is the eigenfunction corresponding to the eigenvalue p 2 . Note that 
the subscript p in (f> p (x) denotes symbolically a set of quantum numbers yet 
to be determined. In the chiral representation in which E 3 and 7 s are both 
diagonal with eigenvalues a — ±1 and % = ±1, respectively, the eigenfunctions 
<f> p (x) have the general form 

(j) p (x) = E pa (x) u ax , (2.4) 

where lo ux are bispinors which are the eigenvectors of E 3 and 7 5 . The scalar 
functions E pa (x) are found to be given by 

E pa (x) = N{n) e l ^ x ° + P^ 2+ P^D n {p), (2.5) 

where N(n) = (4vr|eif|) 1/4 /v / n! is a normalization factor, and D n (p) denotes 
the parabolic cylinder functions [49] with argument p = J2\eH\ {x\ —p2/eH) 



In this article we use the convention in which the metric has the signature g^ u = 
diag(— 1, 1, 1, 1), the Dirac matrices obey {7^,7^} = — 2</ /x " and 7 s = i7°7 1 7 2 7 3 . 
2 Strictly speaking, this maximal set of mutually commuting operators also contains 
the obvious operators — ido, —182 and — ids, which have been omitted here for the 
sake of presentational simplicity. 
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and nonnegative integer index n — 0, 1, 2, . . . given by 

n = / + |sgn(e/J)-i (2.6) 

The nonnegative integer / in Eq. (2.6) labels the Landau levels with / = 
being the lowest Landau level (LLL). Note that in the absence of an external 
electric field, the scalar functions E pa (x) do not depend on the chirality x 
(hence the notation used). 

The eigenvalue p 2 of — (7 • tc) 2 is identified as the momentum squared of a 
massless fermion in an external magnetic field. This can be understood by 
writing (7 • n) 2 = ■ ii\\) 2 + ( r y± ■ 7Tj_) 2 and noticing the following properties: 

-(7|| • 7T||) 2 <f> p {x) = -(7% + 7 3 7T 3 ) 2 M X ) 

= -(vr 2 -nl)(j) p (x) 

= -(p 2 -pI)M x ) 

= pjj0 p (x), (2.7) 

-(7-L • V"_l) 2 (frp(x) = -(7 1 7Ti + 7 2 7T 2 ) 2 (j) p {x) 

= (A + 7T 2 - eHY?) (f) p (x) 
= 2\eH\l(f) p (x) 

= P±M*)> (2-8) 

where, here and henceforth, the subscript || (_L) refers to the longitudinal: 
(j, — 0, 3 (transverse: /i — 1,2) components. Hence, one finds 

2 2,2 
V = P\\ + P± 

= -p 2 +p£ + 2\eH\l, (2.9) 
with / being the quantum number of the quantized transverse momentum 



squared and y2\eH\ the energy gap between adjacent Landau levels (referred 
to henceforth as the Landau energy). 

It is important to note that Eq. (2.6), together with the nonnegativeness of n, 
imposes a constraint on the allowed value of a when / = 0, namely, 

_Jsgn(eiJ) for/ = 0, 

\±sgn{eH) for/>0. [ UJ 

This is one of the subtle points that has been overlooked in the previous 
literature [20,22,23,48] that utilizes the Ritus E p functions. Physically, the 
constraint on a when I = means that the spin of the LLL fermions is always 
aligned parallel to the external magnetic field. 

Following Ritus [46,47], we construct the E p functions as the matrix of the 
simultaneous eigenfunctions (eigenvectors) of the maximal set of mutually 
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commuting operators {(7 • it) 2 , £ 3 ,7 5 }, viz, 

E p (x)= J2' E pa (x)A(a), (2.11) 



a=±l 



where a prime on the summation symbol means that it is subject to the 
constraint (2.10), and 

A(cr) = diag(<5 CT i, 5 al , 5 CT _i) 

= i(l + aS 3 ). (2.12) 

Note that the subscript p in E p (x) denotes symbolically the set of quantum 
numbers {po,P2,P3,l}- 

It is straightforward to verify that (no summation over repeated indices is 
implied) 

A (a) + A(-a) = 1, A(a) A(a') = 5 aa ,A(a), (2.13) 

hence the matrices A(er) can be identified as the projection operators on the 
fermion states with the spin parallel (a = 1) and antiparallel (er = — 1) to the 
external magnetic field. Furthermore, A (a) satisfy the following important 
properties 

[A((7),7j , ] = 0, A(<7)t£ = t£A(-<7), (2.14) 
which are the key relations that will be frequently used in our calculation. 
Using the orthogonal property of the parabolic cylinder functions [49] 

/oo 
dpD n ,{p) D n { P ) = V2^n\ 6 nn ,, (2.15) 
-00 

we find that the E p functions are orthonormal and complete, namely, 

d A xE p (x)E p ,(x) = (27r) 4 5 {4 \p-p')U{l), (2.16) 
f^d A pE p {x)E p {y) = (2n) 4 5^(x-y), (2.17) 

where E p = 7°£^7°. Here we have used the following shorthand notation 

n(0 = (A[ S g„(etf)] for ' = °- (2.18) 
\ 1 for I > 0, V ; 

p 00 poo poo poo 

X^ = E/ d Po dp 2 dp 3 , (2.19) 

J £_q J — CO J —CO J —CO 

5^{p-p>) = S lll 5(p -p' )5(p 2 -p' 2 )5(p 3 -p' 3 ). (2.20) 
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We emphasize that the orthonormal condition of the E p functions for the 
LLL fermions (i.e., I — 0) in Eq. (2.16) differs from that quoted in the lit- 
erature [20,22,23,48], which is applicable only for fermions occupying higher 
Landau levels (i.e., / > 0). The presence of the projection operator A[sgn(eif)] 
in the orthonormal condition of the E p functions for the LLL fermions has an 
important physical consequence that is responsible for an effective dimensional 
reduction in the dynamics of fermion pairing in a strong external magnetic 
field [17,19]. 

Since the E p functions form a complete set of orthonormal Dirac matrix- 
valued functions, they can be used as a basis for the Hilbert space of the 
(perturbative) asymptotic massless fermion states in the presence of a constant 
external magnetic field. Throughout this article we will refer to the space 
spanned by the E p functions simply as the momentum space whenever no 
confusion may arise. 

Explicit calculation shows that the E p functions satisfy an important property 

( 1 -n)E p (x) = E p (x) 1 -p, (2.21) 

where 

P» = (po,0,-sgn(ei7)^2|etf|Z,p 3 ). (2.22) 
For notational simplicity, in Eq. (2.21) we have used the same notation p to 
denote both symbolically the set of quantum numbers {po,P2,P3,l} on the 
left-hand side (LHS) in E p (x), as well as literally the momentum p M on the 
right-hand side (RHS) in 7 • p. However, we emphasize that the quantum 
number p 2 should not be confused with the component of the momentum p 
in the ^-direction, which is determined by the quantum number I as can be 
seen from Eq. (2.22). 

We note that 7 • 7r does not commute with S 3 and 7 5 , hence it is not diagonal 
in the basis spanned by the E p functions, as is evident from Eq. (2.21). 

Upon left multiplying Eq. (2.21) by E p r(x) and integrating over x, one obtains 

J d'xE^^-^E^x) = (27r) i 5^(p-p')U(l) 1 -p. (2.23) 

Eq. (2.23) reveals clearly the advantage of the E p functions in the studies of 
QED in a constant external magnetic field: When expressed in the momentum 
space spanned by the E p functions, the (bare) inverse propagator for fermions 
in a constant external magnetic field is rendered formally identical to that in 
the absence of external fields (up to the projection operator A[sgn(eif)] that 
accounts for the spin alignment of the LLL fermions). 

We are now in a position to show that in the momentum space spanned by 
the E p functions, Eq. (2.1) is rendered formally identical to the Dirac equation 
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for a massless fermion in the absence of external fields. This is achieved by 
considering the Dirac spinor in momentum space ip{p) defined by 

^)=$^E P (*)HP), (2-24) 

where the momentum argument p in ip{p) is given by Eq. (2.22). It follows 
from Eq. (2.21) that in momentum space Eq. (2.1) reads 

( 7 .p)^(p) = 0, (2.25) 

which, as advertised in the beginning of this section, is formally the Dirac 
equation for a massless fermion in the absence of external fields, and whose 
solutions are well known. The condition that Eq. (2.25) has nontrivial solutions 
is given by p 2 = 0, which is exactly the mass shell condition for a massless 
particle. We emphasize that because the E p functions are not solutions of the 
free field Dirac equation (2.1), they do not correspond to the (perturbative) 
asymptotic states of massless fermions in the presence of a constant external 
magnetic field. Nevertheless, as clearly displayed in Eqs. (2.24) and (2.25), 
the latter can be expanded in the basis spanned by the E p functions in terms 
of the usual Dirac (or Weyl) spinors for massless fermions in the absence of 
external fields. 

So far we have been considering only the noninteracting theory for massless 
fermions in a constant external magnetic field. Once the interaction with the 
(quantum) photon field is taken into account, the fermion will receive radia- 
tive corrections, which in turn give rise to the self-energy. Consistent with 
Eq. (2.21), the fermion self-energy in coordinate space E(x, has to satisfy 
a similar property, 

J d A x' x') E p (x') = E p (x) £(p), (2.26) 

where E(p) is a Dirac matrix- valued function (yet to be determined by explicit 
calculation) with the momentum argument p given by Eq. (2.22). 

Following the same steps that lead to Eq. (2.23), one finds the fermion self- 
energy in momentum space to be given by 

E(p,p') = J d A xd A x'E p (x)^(x,x')E pl {x') 

= {2ii) A 5 {A \p-p l )T[{l)Y 1 (p). (2.27) 

With Eqs. (2.23) and (2.27), one finds that the full inverse fermion propagator 



9 



in momentum space is given by [see Eq. (3.1)] 

G-^j/) = J d 4 xd 4 x'E p {x) G-\x,x')E p ,(x') 

= (2vr) 4 S^(p - p>) U(l) [7 • P + E(p)], (2.28) 

which in turn implies that the full fermion propagator in momentum space 
takes the form 

G(p,p') = J d 4 xd 4 x'E p {x)G{x,x')E p ,(x') 

= (2n) 4 S^(p-p')n(l) 1— -. (2.29) 

7 -P + S(p) 

Eqs. (2.27)-(2.29), together with the orthonormal condition Eq. (2.16) of the 
E p functions, imply that in coordinate space one has 

^ = $ E ^ X ) S (P) ^(*0> ( 2 - 3 °) 
G- 1 ^, x>) = ^-0^ 4 E p (x) [7 • p + E(p)] ^(x'), (2.31) 

The set of equations (2.27)-(2.32) constitutes the main ingredients for the 
calculation presented in Sees. 3 and 4. 

We emphasize that the fermion self-energy E(p) in general is noi a diagonal 
matrix, which is contrary to what has been stated in the literature [20,22,23,46- 
48]. We shall now clarify this confusion. It was argued [20,22,46-48] that in 
the presence of a constant magnetic field the fermion self-energy operator £ 
is a Dirac matrix-valued scalar function of the form 

S = S[ 7 -7T,(F^) 2 ], (2.33) 

where = d^A u — d u A fl is the external field strength. In order to highlight 
the fact that the self-energy is in general an off-shell quantity, we have used 
the term self-energy operator instead of the usual term mass operator, which 
can be misleading. 

A comment here is in order. Based on symmetry arguments [46,47], for con- 
stant external fields the self-energy operator E in general depends on two ex- 
tra operators: 7 s F^F^ and a^F^, where F^ v = (1/2) e^ ap F ap and a*" = 
(i/2)[7 M , 7"]. The former vanishes identically in the absence of an external 
electric field. The latter renders either the fermion self-energy or the position 
of the fermion pole (i.e., the fermion dispersion relation) explicitly dependent 
on the sign of the external magnetic field strength, hence is forbidden by a 
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residual Zi symmetry of the rotational 0(3) symmetry that is broken in the 
presence of an external magnetic field (see Sec. 3 for a related discussion). 

Since (F^ti u ) 2 commutes with (7 • 7r) 2 , the self-energy operator E commutes 
with (7 • 7r) 2 and hence is diagonal in the basis spanned by the eigenfunctions 
of (7 • 7r) 2 . While this statement is correct in itself, it does not imply that E 
is diagonal in the basis spanned by the E p functions. This can be understood 
as follows. We note that even though E commutes with (7 • 71) 2 , it does not 
commute with E 3 and 7 5 . Therefore, E is not diagonal in the basis spanned 
by the E p functions, which are constructed in terms of the simultaneous eigen- 
functions (eigenvectors) of the maximal set of mutually commuting operators 
{(7 • 7r) 2 , E 3 , 7 5 }. This is precisely the same reason why the operator 7 • tt, in 
spite of commuting with (7 • 7r) 2 , is not diagonal in the basis spanned by the 
E p functions. 

Last, but not least, we note that E does not anticommute with 7 s . This in 
turn suggests, as one would already have expected, dynamical chiral symme- 
try breaking in the full theory of QED in the presence of a constant external 
magnetic field. While can be confirmed only by explicit nonperturbative calcu- 
lation, the above symmetry arguments seem to indicate that dynamical chiral 
symmetry breaking should be a generic phenomenon independent of the num- 
ber of fermion flavors. Indeed, it is one of the goals of this article to verify in 
a gauge independent manner that in a constant external magnetic field chiral 
symmetry is dynamically broken regardless of the number of fermion flavors. 



3 Schwinger-Dyson equations and the bare vertex approximation 

The Schwinger-Dyson (SD) equations in massless QED in an external mag- 
netic field are well-known in the literature (for reviews see Refs. [50,51]). The 
equations for the full fermion propagator G(x,y) are given by 3 

G- 1 (x,y) = S- 1 (x,y) + J:(x,y), (3.1) 
Z(x,y)=te 2 J ^x'd'y'YG^x'^ix'^V^y'), (3.2) 

where S(x, y) is the bare fermion propagator in the external gauge field 
given by Eq. (2.2), E(x, y) is the fermion self-energy and T u (x, y, z) is the full 



3 We note that in writing the fermion self-energy S(ar, y) as in Eq. (3.2), it is 
implicitly assumed that the vacuum current vanishes. In Appendix A we show that 
the vacuum current does indeed vanish identically in a constant external magnetic 
field. 
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i/WVWVr = \A/WWV + 




Fig. 1. The SD equations in QED in an external magnetic field. The thin (heavy) 
line with an arrow denotes the bare (full) fermion propagator in the presence of an 
external magnetic field. The thin (heavy) wiggly line denotes the free (full) photon 
propagator. The vertex with(out) a heavy dot denotes the full (bare) fermion-photon 
interaction vertex. 

vertex. The full photon propagator V^ u (x,y) satisfies the equations 



where D^ u (x,y) is the free photon propagator (defined in covariant gauges) 
and U fiU (x, y) is the vacuum polarization. The diagrammatical representation 
of the SD equations (3.1)- (3. 4) is depicted in Fig. 1. 

Since the dynamics of fermion pairing in a strong magnetic field is dominated 
by the lowest Landau level (LLL) [18-22,24,25], we will consider the prop- 
agation of, as well as radiative corrections originating only from, fermions 
occupying the LLL. This is referred to as the lowest Landau level approxi- 
mation (LLLA) in the literature [18-29]. Therefore, except where otherwise 
explicitly stated, the fermion propagator and the fermion self-energy in the 
SD equations (3.1)- (3.4) as well as in the rest of this article will be taken to 
be those for the LLL fermions. For notational simplicity, no separate notation 
will be introduced. 

It is well known that the SD equations (3.1)- (3.4) correspond to an infinite 
hierarchy of integral equations. This is because the full vertex (i.e., the three- 
point vertex function) depends on four-point vertex functions which in 
turn depend on higher-point vertex functions, and so on. In order to reduce 
the SD equations to a closed system of integral equations that is tractable, a 
truncation scheme of the SD equations needs to be employed by truncating 
this infinite hierarchy at some point. Diagrammatically, a truncation scheme 
of the SD equations corresponds to a resummation of a selected infinite subset 
of diagrams arising from every order in the loop expansion. The simplest 
truncation is done at the level of the three-point vertex function by expressing 




(3.3) 
(3.4) 



ie 2 tr J d A x'd 4 y' 7a* G(x, x')T v (x', y', y) G(y f , x) 
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Fig. 2. The SD equations in the bare vertex approximation. 

the full vertex T M in terms of the (bare or full) entities that already appeared 
in the SD equations. 

To this end we will work in the bare vertex approximation (BVA), in which 
the vertex corrections are completely ignored. This is achieved by replacing 
the full vertex in the SD equations (3.1)-(3.4) by the bare one, viz, 

r»(x,y,z) =^5 {i \x-z)5 {i \y-z). (3.5) 

The diagrammatical representation of the resultant SD equations truncated 
in the BVA is depicted in Fig. 2. This truncation is also known as the (im- 
proved) rainbow approximation and has been employed extensively in the lit- 
erature [18-20,22,24-29]. However, we emphasize that unlike what has usually 
been done in the literature, here we will not confine ourselves to a particular 
gauge (usually the Feynman gauge) [18-20,22], nor will we make the assump- 
tion that the BVA (3.5) is valid only in a certain gauge [24-29]. Instead, it 
is our aim to prove that the BVA (3.5) is a consistent truncation of the SD 
equations (3.1)- (3.4) within the LLLA. The dynamical fermion mass, obtained 
as the solution of the truncated SD equations evaluated on the fermion mass 
shell, is manifestly gauge independent. In the weak coupling regime that we 
consider, such a gauge independent approach allows one to resum consistently 
an infinite subset of diagrams that arise from every order in the loop expan- 
sion and whose contributions are of leading order in the gauge coupling, thus 
leading to a consistent and reliable calculation of the dynamical fermion mass. 

The main ingredient in the proof of the gauge independence of physical quan- 
tities is the Ward-Takahashi (WT) identity satisfied by the vertex and the in- 
verse fermion propagator. It is known [13,52] that in ordinary QED the BVA 
is not a consistent truncation of the SD equations because the correspond- 
ing WT identity between the full fermion propagator and the bare vertex in 
general cannot be fulfilled. The situation however changes drastically in the 
presence of a strong, constant external magnetic field as the motion of low- 
energy fermions is restricted in directions perpendicular to the magnetic field, 
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leading to an effective dimensional reduction from (3 + 1) to (1 + 1) in the 
dynamics of fermion pairing in a strong magnetic field. 

The WT identity in the BVA within the LLLA was first studied in Ref. [23]. 
It was shown that in order to satisfy the WT identity in the BVA within the 
LLLA, the LLL fermion self-energy in momentum space has to be a momen- 
tum independent constant. However, due to a mistake in Ref. [23] regarding 
the Dirac matrix structure in the orthonormal condition of the Ritus E p func- 
tions for the LLL fermions [see Eq. (2.16) and the discussion that follows], the 
calculation therein requires further investigations. As we will see below, cal- 
culation performed by using the correct orthonormal condition shows that the 
conclusion obtained in Ref. [23] on the WT identity in the BVA remains valid. 
The reliability of such a momentum independent approximation (in that the 
fermion self-energy has very weak momentum dependence) and, consequently, 
of the WT identity in the BVA has been verified within the LLLA in a certain 
gauge in the momentum region relevant to the dynamics of fermion pairing in 
a strong magnetic field [24-28] (see further discussions in Sec. 5). 

We now derive the WT identity in the BVA both within and beyond the LLLA. 
This allows us to (i) correct a few oversights found in Ref. [23] , (ii) verify the 
conclusion obtained therein on the WT identity in the BVA within the LLLA, 
and (iii) show that the WT identity in the BVA can be satisfied only within 
the LLLA. It is noted that in the rest of this section the fermion propagator 
and the fermion self-energy will in general not be restricted to those for the 
LLL fermions. 

Following the standard procedure for deriving the WT in QED [53,54,23], one 
finds the WT identity in the BVA in the presence of an external gauge field 
to be given by 

-d^[5 {4 \x-z)5 {4 \y-z)Y} =t[5 {4) (y-z) - 8^(x - z)]G-\x,y), (3.6) 

where it is noted that G(x, y) is the full fermion propagator in an external 
gauge field. Since the inverse fermion propagator in the above equation does 
not depend on z, performing the Fourier transform on the variable z leads to 

e -iq-x S W( X _ y ) 1 . q= ( e -i9-* _ e -w-w) G~ 1 (x, y), (3.7) 

where is the momentum carried by the photon. Upon left and right multi- 
plying Eq. (3.7) by E p (x) and E p '(y), respectively, and integrating over x and 
y, we obtain the following identity 

J d 4 xe-^ x E p (x) r qE p ,(x) = J d A x d A y (e~ iq ' x -e~ iq ' y ) E p (x) G~ l (x, y) E pl (y), 

(3.8) 

where the subscripts p and p' denote the sets {po,P2,P3,l} and {p' , p' 2 , p' 3 , 1'} , 
respectively. 
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Let us first derive the WT identity in the BVA within the LLLA, namely, 
I — I' — in Eq. (3.8). Using the properties of the E p functions [see Eqs. (2.11), 
(2.12) and (2.5)], one finds that the LHS of the WT identity (3.8) is simplified 
to 

J d A xe- iq - x E p (x)-i-qE p ,(x) = (2tt) 4 5 {3 \p + q - p') e^Wl e -m( P 2+ P ' 2 )/^H 

x A[sgn(e#)] 7 || • q\\ A[sgn(e#)], 

(3.9) 

where the transverse components 7j_ • q± decouple due to the second property 
of A(a) in Eq. (2.14). Here q\ = q\ + q\ and 

<S(3)(p + q _ p ') = 5 ( po + go _ p/ o ) 5 ( p2 + g2 _ p2 ) 5 ( p3 + gs _ p ' 3 y (3.io) 

Likewise, the RHS of the WT identity (3.8) can be simplified to 

J d 4 x d 4 y (e-^ x - e-**) E p {x) G~\x, y) E p ,(y) 

= (2vr) 4 5 (3) (p + q- p') e-^ eH \ e -m( P 2 +P ' 2 )/2eH 

x A[sgn(eiJ)] [7,, • (p> - p),, + E(p\) - S(p M )] A[sgn(e//)], (3.11) 

where used has been made of Eq. (2.31), and the momenta p and p' on the RHS 
in the square brackets are given by p^ = (p ,0,0,p 3 ) and p'^ = (p' ,0,0,p' 3 ), 
respectively [see the remarks after Eq. (2.22)]. Upon collecting the above re- 
sults and using the delta functions to express p'^ in terms of (p + q)\\, we can 
rewrite the WT identity (3.8) in the following compact form 

A[sgn(e#)] {7||-?||-[7ir(P+9)|| + S(P||+9||)] + [7irP|| + S(p||)]} A[sgn(e/f)] = 0. 

(3.12) 

In the above expression, we have purposely kept the tree-level terms, which 
cancel among themselves, to highlight the most general form of the WT iden- 
tity in the BVA within the LLLA. 

To proceed further, one needs to know the Dirac matrix structure of the 
fermion self-energy Based on symmetry arguments [23], for a constant 

external magnetic field E(p) takes the following form 

E(p) = A\\(p) 7|| • P\\ + A ± (p) 7 _l • p ± + B(p), (3.13) 

where A»±{p) and B(p) are functions of the longitudinal and transverse mo- 
mentum squared as well as the magnitude of the external magnetic field, yet 
to be determined by explicit calculation. 

A comment here is in order. It was shown in Ref. [23] that for a constant 
external magnetic field along the ^-direction, the most general structure of 
contains two extra terms that are proportional to if E 3 , where H is the 
strength of the constant external magnetic field. Detailed analysis indicates 
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that such terms, if present, lead to either the fermion self-energy or the position 
of the fermion pole being explicitly dependent on the sign of the external 
magnetic field strength H, hence are forbidden by a residual Z 2 symmetry 
(i.e., H — > —H) of the rotational 0(3) symmetry that is broken in the presence 
of an external magnetic field (see Sec. 2 for a related discussion). 

Upon substituting Eq. (3.13) into Eq. (3.12), one finds that the transverse com- 
ponents proportional to 7j_ • p± in Eq. (3.13) decouple from the WT identity. 
Furthermore, the requirement that the WT identity is fulfilled for arbitrary 
momenta p and q regardless of the sign of H, leads to the conclusion that 
S(py) has to be a momentum independent constant. 

Before ending this section, let us derive the WT identity in the BVA but beyond 
the LLLA, where one of the fermions occupies the lowest Landau level while 
the other fermion occupies a higher Landau level. In other words, we consider 
Eq. (3.8) with I = and V > 0. We shall show that the WT identity in the BVA 
can be satisfied only within the LLLA. This is an important conclusion that 
allows for an unequivocal refutation of the results and conclusions obtained 
in Ref. [29], which can be attributed to an inconsistent truncation of the SD 
equations beyond the LLLA (see Sec. 5 for details). 

Following similar steps as above, one can rewrite the WT identity (3.8) in the 
BVA but beyond the LLLA in the following compact form (V > 0) 

A[sgn(e#)] {K v (ql, <p) ^ • q {l - 7 • p' + 7, • p\\ - E(p') + E(p||) 

+ K l ,_ 1 (ql,<p)>y ± .q ± } = 0, (3.14) 

where use has been made of the properties of A(cr) in Eq. (2.14) as well as 
the Dirac matrix structure of the fermion self-energy in Eq. (3.13), and the 
momenta p, p' and q are respectively given by [see the remarks after Eq. (2.22)] 

= (Po,0,0,p 3 ), 



pj, = W A-8g*(eH)y/2\eH\l',j/ 3 ), (3.15) 
% = (p'o -Po,qi,p' 2 -P2,p' 3 ~Ps)- 



In the above expression, we have use the shorthand notation 

1/2] v 

(3.16) 



Kv(q 2 ± ,ifi) = -j= exp[-i sgn(eif) 1'ip] 



2 \ V2l 



iSgn{eH \^H\ 

where f = arctan(g 2 /o , i)- It is an easy exercise to verify that for either 
a momentum independent or a momentum dependent fermion self-energy, 
Eq. (3.14) cannot be satisfied for arbitrary momenta p and p' regardless of the 
sign of H. This result is not unexpected because, as is evident from Eq. (3.14), 
there is no longer an effective dimensional reduction in the dynamics of fermion 
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pairing beyond the LLLA. The situation is then similar to that in ordinary 
QED, in which it is known that the BVA is not a consistent truncation of the 
SD equations [13,52]. Consequently, we conclude that (i) the WT identity in 
the BVA can be satisfied only within the LLLA with a momentum indepen- 
dent fermion self-energy; (ii) in order to go beyond the LLLA one has to use 
truncation schemes of the SD equations that consistently account for vertex 
corrections. To the best of our knowledge, such a consistent truncation of the 
SD equations beyond the LLLA has not appeared in the literature. 



4 On-shell gauge independence in the bare vertex approximation 

As per the WT identity in the BVA and within the LLLA (3.12), we can write 
the self-energy for the LLL fermion as 

S(p||)=m(0, (4.1) 

where m(£) is a momentum independent but gauge dependent constant, with £ 
being the gauge fixing parameter in covariant gauges. It is important to note 
that m(£) depends implicitly on £ through the full photon propagator 
in Eq. (3.2). We emphasize that because of its implicit ^-dependence, m(£) 
should not be taken for granted to be the dynamical fermion mass, which is 
a gauge independent physical quantity. 

We now begin the proof that the BVA is a consistent truncation of the SD 
equations (3.1)-(3.4), in which m(£) is ^-independent and hence can be iden- 
tified unambiguously as the dynamical fermion mass, if and only if the trun- 
cated SD equation for the fermion self-energy is evaluated on the fermion mass 
shell. 

We first recall that, as proved in Ref. [15], in gauge theories the singularity 
structures (i.e., the positions of poles and branch singularities) of gauge boson 
and fermion propagators are gauge independent when all contributions of a 
given order of a systematic expansion scheme are accounted for. Consequently, 
this means the dynamical fermion mass has to be determined by the pole of 
the full fermion propagator obtained in a consistent truncation scheme. 

In momentum space the full propagator for the LLL fermion is given by [see 
Eq. (2.29)] 

G(p\\) = A[sgn(etf)], (4.2) 

7|| -P\\ + S(P||) 

where, as noted in Sec. 2, A[sgn(eiJ)] is the projection operator on the fermion 
states with the spin parallel to the external magnetic field. Assume for the 
moment that the BVA is a consistent truncation of the SD equations in the 
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LLLA, such that the position of the pole of the LLL fermion propagator G(p») 
is gauge independent. It follows that in accordance with the WT identity in the 
BVA (3.12) or, equivalently, the momentum independence of the LLL fermion 
self-energy, we have 

£(pil) = £(pjj = -m 2 )=m, (4.3) 

where m is the gauge independent, physical dynamical fermion mass, yet to 
be determined by solving the truncated SD equations self-consistently. With 
the LLL fermion self-energy given by Eq. (4.3), the WT identity in the BVA 
(3.12) reduces to a tree-like form 

7|| • q\\ = (7|| • P\\ +m)- [7|| • (p - g)|| + m\. (4.4) 



What remains to be verified in our proof that the BVA is a consistent trun- 
cation of the SD equations in the LLLA is the following statements: (i) the 
truncated vacuum polarization is transverse, (ii) the truncated fermion self- 
energy when evaluated on the fermion mass shell (p| = — m 2 ) is manifestly 
gauge independent. We highlight that the fermion mass shell condition is one 
of the most important points that has gone unnoticed in the literature, where 
the truncated fermion self-energy used to be evaluated off the fermion mass 
shell at, say, p\ = [18-20,22,24-29]. 

Let us first consider the vacuum polarization, which in the BVA reads 

IF* (x, y) = -te 2 tr ^ G(x, y)^G{y,x), (4.5) 

where the LLL fermion propagator G(x,y) is given by Eq. (2.32) with the 
summation over / restricted to the I = term and, as per the WT identity in 
the BVA (3.12), a momentum independent self-energy given by Eq. (4.3). In 
other words, one has 

where a prime on the summation symbol means the summation over / is re- 
stricted to the / = term. Taking the Fourier transform of Eq. (4.5) and using 
the properties in Eq. (2.14), we find the vacuum polarization in momentum 
space to be given by 

n-(g) = --|etf|exp(-^U tr/^L 
W 2tt 1 1 P \ 2\eH\J J (2tt) 2 

x 7i| / \ ■ A[sgn( e g)]. (4.7) 

11 7|| • p\\ + m 11 7|| • (p — q)\\ + m 

The presence of A[sgn(eif)] in Eq. (4.7) is a consequence of the LLLA, which, 
as explicitly displayed in Eq. (4.7), leads to an effective dimensional reduction 
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from (3 + 1) to (1 + 1) as the LLL fermions couple only to the longitudinal 
components of the photon field. 

The WT identity in the BVA (4.4) guarantees that the vacuum polarization 
n^(g) is transverse, viz, 

g"lWg) = 0. (4.8) 

Explicit calculation in dimensional regularization shows that the 1/e pole cor- 
responding to an ultraviolet logarithmic divergence cancels, leading to 

U^(q) = Tl(^)(gf-^f]. (4.9) 

This in turn implies that the full photon propagator takes the following form 
in covariant gauges (£ = 1 is the Feynman gauge): 

The polarization function n(gjj, g 2 ) is given by 

where a = e 2 /47r is the fine-structure constant and 

i x /i + i/u + 1 

F « = i T—r lo s \ = — • ( 4 - 12 ) 

2Uyj\ + I/U yfl + I/U ~ 1 

Our result for n(g|, g 2 ) agrees with those obtained in the literature [19,21,55]. 
It is worth noticing that the result in Ref. [21] is valid only in the limit m — > 0. 

A detailed analysis of the analytic structure of the function F(u) shows that 
F{u) has an imaginary part for u < —1, viz, 

lmF(u) = n Q(-u-l), (4.13) 

2Uy/l + I/U 

where 6(x) is the Heaviside step function. Hence the polarization function 
il(gjj, g 2 ) is complex for gjj < —4m 2 , with an imaginary part given by 

Imn(g 2 ,gi) = \eH\ sgn(g ) exp(-^L) hnF^-^j. (4.14) 

We note that the extra factor of — sgn(g ) is due to the analytic continuation 
to the complex g°-plane, namely, g° — > q° + ie with e — > + . As will be 
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Fig. 3. Plot of the real and imaginary parts of F(u) as a function of u. 

seen momentarily, the property that lmll(qpq\) is an odd function of go 
has an important consequence on the dynamical fermion mass. The real and 
imaginary parts of F{u) are depicted in Fig. 3. The real part of F(u) has the 
following asymptotic behavior: 



ReF(u) 



f + 0(u 2 ) 



1 

2u 



log(4H) + C?| 



for \u\ <C 1, 
for \u\ ^> 1. 



(4.15) 



The polarization effects modify drastically the propagation of virtual photons 
in a constant external magnetic field. While photons of momenta <C m 2 
remain unscreened, photons of momenta m 2 <C \q 2 \ <C \eH\ and q\ <C \eH\ 
are screened with a characteristic length L = (2a\eH\/n)' 1 / 2 . It is noted 
that the upper limit \eH\ in the range of is due to the LLLA. As will be 
seen in Sec. 5, this screening effect renders the rainbow approximation [18- 
20,22], in which the bare vertex as well as the free photon propagator are used, 
completely unreliable in this problem. 



We next consider the fermion self-energy, which in the BVA reads 
E(x, y) = ie 2 ^ G(x, y) Y T>^(x, y), 



(4.16) 



where G(x, y) is the LLL fermion propagator given by Eq. (4.6) and V fMU (x, y) 
is the full photon propagator 



(2k) 



(4.17) 



with V^(q) given by Eq. (4.10). 
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Following similar steps in Sec. 3 in the derivation of the WT identity by left and 
right multiplying Eq. (4.16), respectively, by E p (x) and E p >(y) corresponding 
to the LLL fermion (i.e., I — I' — 0), and integrating over x and y, we find 
that in momentum space the fermion self-energy evaluated on the fermion 
mass shell is given by 



m A[sga(eH)] —ie z 



d 4 q 



exp 



2\eH\ J " 7|| • (p — g)|| + m 



7y 



xV^(q) A[sgn(eH)} 



(4.18) 



p\=-m? 



where use has been made of Eq. (2.14). The presence of A[sgn(eiJ)] in (4.18) 
is again a consequence of the LLLA. 



We note that the RHS of Eq. (4.18) appears to contain a gauge dependent 
part that arises from the gauge dependent term in T>^ u {q) [see Eq. (4.10)]. 
The WT identity in the BVA (4.4) guarantees that this would-be gauge de- 
pendent contribution to the fermion self-energy (denoted symbolically as 
is proportional to ■ p\\ +m). An explicit calculation shows that is given 
by 



E € = a - 1) (7|| • P|| + m) J q dx(l-x)J-^± exp [-^ 



x (1 ~ x)qj - xm( lv Pll - m) 

X [(l- x )ql + x(l-x)pl + xm*]> A ^ n ^^' ^ iy ^ 



Since is proportional to -p\\ +m), it vanishes identically on the fermion 
mass shell p| = — m 2 or, equivalently, 7y ■ p\\ +m — 0. This, together with the 
transversality of the vacuum polarization, completes our proof that the BVA 
is a consistent truncation of the SD equations. Consequently, the dynamical 
fermion mass, obtained as the solution of the truncated SD equations evaluated 
on the fermion mass shell, is manifestly gauge independent. 



Having proved the gauge independence of the dynamical fermion mass in the 
BVA, we are now ready to find m by solving Eq. (4.18) self-consistently. Note 
that the transverse components in T>^ v (q) decouple in the LLLA. Following 
the same argument given above in the proof of the on-shell gauge indepen- 
dence of the fermion self-energy, one can verify fairly easily that contributions 
from the longitudinal components in T>^ u (q) proportional to q^q\\/q 2 are also 
proportional to (711 -pn +m), hence they vanish identically on the fermion mass 
shell. Therefore, only the first term in T>^ v (q) proportional to contributes 
to the on-shell SD equation (4.18). Consequently, the Dirac matrix structures 
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on both sides of Eq. (4.18) are consistent. With these remarks we obtain 

o f d 4 q 2m ( q\ \ 

m = -it \ exp 1 ' ' 



x 



(27r) 4 (p-g)J + m 2 L \ 2\eH 
q 2 + ReU(q 2 q 2 ± )-tlmU(q 2 ql) 



[g2 + Ren(g2 g 2 )]2 + [Imn(g 2 



(4.20) 



Since ImIT(gy, qj_) is an odd function of go [ see Eq. (4.14)], the imaginary part 
of the above integral vanishes identically. As a consequence, the dynamical 
fermion mass m determined by the on-shell SD equation (4.20) is a real- valued 
quantity as it should be. Anticipating that the dynamical fermion mass is much 
less than the Landau energy, i.e., m 2 <C \eH\, we find that the real part of the 
above integral is dominated by the region of momentum m 2 <C \q 2 \ <C \eH\ 
and q\ <C \eH\. Therefore, the imaginary part of the polarization function 
ImiT(g|, q\) in the denominator in Eq. (4.20) can be neglected. Using the 
mass shell condition py 1 = (m, 0) that corresponds to a LLL fermion at rest 
and performing a Wick rotation to Euclidean space, we find that Eq. (4.20) 
in Euclidean space reads 

a r 2 m f°° 2 exp(—qj_/2\eH\) 

m= 2^J dqil q 2 + (q 4 -my + m2Jo Q± qf + q\ + ReU(q 2 , q \) ' (421) 

where q 2 = q 2 + q\. 

Before proceeding further we note that, like what has usually been done in 
the literature, Eq. (4.21) can be obtained directly from Eq. (4.18) by making 
a Wick rotation to Euclidean space. While, as remarked in Ref. [53], this pro- 
cedure is easy to justify in perturbation theory by neglecting the possibility 
of dynamically generated singularities in the first and third quadrants of the 
complex energy plane, it is highly nontrivial in nonperturbative studies and 
must be performed with care. Our detailed analysis of the photon polarization 
in the BVA shows clearly that no such singularities will be generated dynam- 
ically in the first and third quadrants of the complex g -plane, thus leading to 
a justification for the Wick rotation to Euclidean space in this problem. 

The generalization of our result to the case of QED with Nf fermion flavors 
can be done straightforwardly by the replacement Il{q 2 ) q 2 i) — > NfU(q 2 ,q\) 
in Eqs. (4.10), (4.20) and (4.21). We have numerically solved Eq. (4.21) to 
obtain m as a function of a for several values of Nf. The results are displayed 
in Fig. 4. Numerical analysis shows that the solution of Eq. (4.21) can be fit 
by the following analytic expression: 



m = a J2\eH\ (3(a) exp 



7T 



alog(b/Nf a)\ 

where a is a constant of order one, b ~ 2.3, and (3(a) ~ Nf a. 



(4.22) 
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Fig. 4. Plot of m as a function of a for several values of Nf. 



It is instructive to consider the case of QED with a large number of fermion 
flavors, i.e., the 1/Nf expansion in large- Nf QED [21,28]. The advantage of 
the 1/Nf expansion is that it allows for a calculation of the dynamical fermion 
mass without a restriction to small a. Specifically, we take the limit Nf — > oo 
with Nf a finite and fixed. This can be achieved by a rescaling of the coupling 
constant e — > e/^Nf, which together with Eq. (4.22) leads to the following 
result for the dynamical fermion mass in the large- Nf limit: 



where a is the rescaled fine-structure constant, which remains finite and fixed 
in the limit Nf — > oo as noted above. We will argue in the next section that 
Eq. (4.23) can also be obtained from the 2PI (two-particle-irreducible) effective 
action truncated at next-to-leading order in the 1/Nf expansion, which resums 
exactly the same infinite subset of diagrams as the SD equations truncated in 
the BVA. 

We end this section by showing that the on-shell gauge independence of the 
fermion self-energy in the BVA can also be understood in terms of the wave 
function renormalization, i.e., the residue at the fermion pole. Since the gauge 
dependent contribution as well as contributions from the terms proportional 
to q^q\\/qf\ in the full photon propagator are proportional to • p\\ + m), 
they contribute only to the wave function renormalization in the full fermion 
propagator if the fermion self-energy were evaluated off the fermion mass shell. 
Hence they do not affect the position of the fermion pole, which is at py = — m 2 
and gauge independent. This is precisely the reason that the physical dy- 
namical fermion mass has to be determined by the pole of the full fermion 
propagator in a consistent truncation. The wave function renormalization is 




(4.23) 
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gauge dependent, but this is perfectly fine since it is not a physical observable. 
Furthermore, the wave function renormalization can be absorbed into a redef- 
inition of the fermion field. This in turn implies that up to a field redefinition 
the physical fermion pole at p 2 = — m 2 has unit residue. 



5 Comparison to previous works and discussions 

Our results for the dynamical fermion mass differ from those obtained in the 
BVA in Refs. [18-20,22,24-29], in which the truncated fermion self-energy were 
evaluated exclusively off the fermion mass shell at p| = 0. Consequently, the 
results for the dynamical fermion mass obtained in these earlier studies are 
inevitably gauge dependent and cannot be identified unambiguously as the 
physical dynamical fermion mass. We now argue that to a large extent those 
earlier results can be attributed to inconsistent truncation schemes, gauge 
dependent artifacts, or both. 

In Refs. [18-20,22] chiral symmetry breaking in a strong magnetic field was 
first studied in the so-called rainbow approximation, in which the bare vertex 
and the free photon propagator were used. Apart from the gauge dependence, 
these earlier results are found to have a very different functional dependence 
on the gauge coupling. This is a consequence of an inconsistent truncation of 
the SD equations in that, as displayed in Fig. 2, the photon propagator that 
enters the SD equations in the BVA is the full one, hence could be replaced 
by the bare one if and only if the vacuum polarization effects are shown to be 
negligible. However, this is not the case for the problem under consideration. 
As discussed in Sec. 4, there is a strong screening effect in the longitudinal 
components of the photon propagator for photons with momenta m 2 <C \q 2 \ -C 
\eH\ and q\ <C \eH\. This is precisely the reason that the numerical values of 
the dynamical fermion mass found in Refs. [18-20,22] tend to be overestimated 
by several orders of magnitude when compared to what we have obtained in 
Sec. 4. 

More recently, chiral symmetry breaking in a strong magnetic field has been 
studied in the so-called improved rainbow approximation [24,25,28,29], in 
which the bare vertex and the full photon propagator were used. We note 
that the improved rainbow approximation used in Refs. [24,25,28,29] is ex- 
actly the same as the bare vertex approximation used in this article. It can be 
verified fairly easily that the truncated SD equations in both approximations 
resum identically the same infinite subset of diagrams. 

The authors of Refs. [24,25,28] claimed that (i) in covariant gauges there are 
one-loop vertex corrections arising from the term q^q\\/q 2 q 2 in the full photon 
propagator that are not suppressed by powers of a (up to logarithms) and 
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hence need to be accounted for; (ii) there exists a noncovariant and nonlocal 
gauge in which, and only in which, the BVA is a reliable truncation of the 
SD equations that consistently resums these one-loop vertex corrections. The 
gauge independent analysis in the BVA, as presented in this article, shows 
clearly that such contributions vanish identically on the fermion mass shell. 
Furthermore, it is evident that if one expands diagrammatically the SD equa- 
tions in the BVA, one finds by induction that to all orders in the loop expansion 
there is not a single diagram with vertex corrections being resummed by the 
SD equations (see Fig. 2). This statement is true in any gauge (be it covariant 
or otherwise), because gauge fixing affects only the explicit form of the (free 
or full) photon propagator but not the topology of the Feynman diagrams. 
Therefore, both the above quoted conclusions are incorrect. 

We emphasize that the WT identity is a necessary condition for establishing 
the gauge independence of the dynamical fermion mass, but it is far from suf- 
ficient. While the WT identity guarantees the truncated vacuum polarization 
is transverse, it guarantees only the truncated on-shell fermion self-energy is 
gauge independent. Hence, the dynamical fermion mass is gauge independent 
only when determined by the position of the fermion pole obtained in a con- 
sistent truncation. This is tantamount to evaluating the truncated fermion 
self-energy on the fermion mass shell. Even though the WT identity in the 
BVA is verified in Refs. [24,25] in a particular noncovariant and nonlocal 
gauge, this does not guarantee that the dynamical fermion mass obtained 
therein from the truncated fermion self-energy evaluated off the fermion mass 
shell will be gauge independent. In fact, the authors of Refs. [24,25,28] simply 
chose a particular noncovariant and nonlocal gauge such that the gauge de- 
pendent contribution happens to cancel contributions from terms proportional 
to q^q\\/qf\ in the full photon propagator [56]. Our gauge independent analysis 
in the BVA reveals clearly that such a gauge fixing not only is ad hoc and 
unnecessary, but also leaves the issue of gauge independence unaddressed. 

In Ref. [29] its authors claimed that (i) in QED with Nf fermion flavors a 
critical number N cr exists for any value of a, such that chiral symmetry re- 
mains unbroken for Nf > N cr ; (ii) the dynamical fermion mass is generated 
with a double splitting for Nf < N cr . As can be gleaned clearly from Fig. 4 
or from the large- Nf result, Eq. (4.23), both the conclusions quoted above are 
incorrect. They are gauge dependent artifacts of an inconsistent truncation: 
On the one hand, the SD equation for the momentum independent fermion 
self-energy was obtained (in an unspecified "appropriate" gauge) in the BVA 
within the LLLA. On the other hand, the vacuum polarization was calculated 
in the BVA but beyond the LLLA. This however is not a consistent truncation 
of the SD equations since, as shown in the end of Sec. 3, the WT identity in 
the BVA can be satisfied only within the LLLA. 

The results of Ref. [29] suggest that in the inconsistent truncation as well 
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as in the unspecified "appropriate" gauge used there, the unphysical, gauge 
dependent contributions from higher Landau levels are so large that they be- 
come dominant over the physical, gauge independent contribution from the 
LLL and lead to the authors' incorrect conclusions. Hence, we emphasize that 
the LLL dominance in a strong magnetic field should be understood in the 
context of consistent truncation schemes as follows. Contributions to the dy- 
namical fermion mass from higher Landau levels that are obtained in a (yet to 
be determined) consistent truncation of the SD equations are subleading when 
compared to that from the LLL obtained in the consistent BVA truncation. 
Our analysis of the WT identity in the BVA but beyond the LLLA indicates 
that in order to go beyond the LLLA, one has to use truncation schemes that 
consistently account for vertex corrections. To the best of our knowledge, con- 
sistent truncations of the SD equations that include vertex corrections either 
within or beyond the LLLA have not appeared in the literature. In view of 
this, a detailed study of the general structure as well as the analytic properties 
of the full vertex akin to that in Ref. [57] but in the presence of a constant 
external magnetic field is of crucial importance. 

Extension to a momentum dependent fermion self-energy in the BVA has 
been studied in a certain gauge by neglecting corrections to the vacuum po- 
larization [24-28]. It was found that the fermion self-energy has very weak 
momentum dependence in the region relevant to fermion pairing in a strong 
magnetic field. Since a momentum dependent fermion self-energy cannot ful- 
fill the WT identity in the BVA, the results obtained in these studies are 
inevitably gauge dependent. It is, however, important to note that when com- 
pared to the corresponding results obtained using a momentum independent 
fermion self-energy, these results lend support to the reliability of a momen- 
tum independent fermion self-energy and, consequently, of the WT identity 
in the BVA in the momentum region <C \eH\ that is relevant to the dy- 
namics of fermion pairing in a strong magnetic field. Similar study has also 
been carried out in a certain gauge in the case of QED with an extra charge 
neutral, self-interacting scalar field that couples to the fermion through the 
Yukawa interaction [58,59]. 

Despite obtained from gauge dependent calculations in a certain gauge, as 
will be argued momentarily, the conclusion from these studies [24-28] that a 
momentum independent fermion self-energy or, equivalently, the WT identity 
in the BVA is reliable in the momentum region relevant to chiral symmetry 
breaking in a strong magnetic field remains valid. In fact, this is precisely the 
underlying physical reason that the BVA is a consistent truncation of the SD 
equations in the LLLA. 

In the BVA the presence of momentum dependent terms in the fermion self- 
energy leads to violation of the WT identity, hence introduces gauge dependent 
contributions to the SD equations even though the latter are evaluated on the 
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fermion mass shell. We now argue that such gauge dependent contributions are 
of higher orders when compared to the dynamical fermion mass obtained from 
a momentum independent fermion self-energy. Therefore, within the LLLA 
the dynamical fermion mass can be reliably calculated in the BVA by using a 
momentum independent fermion self-energy. 

To begin with, we first note that the infinite subset of diagrams being re- 
summed in Refs. [24-28] are identical to, as well as in the same kinematic 
region as, those in the present article. The only difference is that in those 
references the calculation is carried out in a certain gauge and the fermion 
self-energy is evaluated off the fermion mass shell at p| = 0. Since the dy- 
namical fermion mass is much less than the Landau energy, i.e., m 2 <C \eH\, 
the LLL fermion self-energy with momentum <C \eH\ can be ex- 

panded in powers of p 2 /\eH\ around the fermion mass shell in the BVA (i.e., 

A\\(p\\) = a{i)x + 0(x 2 ), (5.1) 
B{p\\) = m[l + b{i)x + 0{x 2 )l (5.2) 

where x — {p\ + m 2 )/\eH\, and a(£) and 6(£) are ^-dependent coefficients 
presumably of order one. This, together with the gauge dependent results 
obtained in the LLLA [24-28] regarding the reliability of a momentum inde- 
pendent fermion self-energy in the BVA in the momentum region \p?,\ <C \eH\, 
indicates clearly that the gauge dependent contributions to the dynamical 
fermion mass arising from the momentum dependent terms in ^4||(p||) and 
B(p\\) are of order |p||/|eif|. Hence physical quantities sensitive only to the 
kinematic region of momenta <C \eH\ can be reliably calculated in the 
BVA, and are gauge independent up to corrections of order |p||/|eif | <C 1. 
This in turn means that, consistent with the underlying assumption of the 
LLLA, the kinematic regime of low-energy photon exchange qj_ <C \eH\) 
is precisely that which determines the dynamics of fermion pairing in a strong 
magnetic field. 

We are now in a position to establish direct contact between the SD equations 
truncated in the BVA and the 2PI (two-particle-irreducible) effective action 
truncated at the lowest nontrivial order in the loop expansion as well as in the 
1/Nf expansion. Since the 2PI effective action is constructed in terms of the 
full propagators and bare vertices, the corresponding WT identities usually 
cannot be satisfied. In general, any truncation of the 2PI effective action is 
inevitably gauge dependent and hence, strictly speaking, an inconsistent one. 
Nevertheless, it has been shown recently in Refs. [60-62] that the truncated on- 
shell 2PI effective action has a controlled gauge dependence with the explicit 
gauge dependent terms always appearing at higher order. This is reminiscent 
of the situation that the dynamical fermion mass obtained from the on-shell 
SD equations truncated in the BVA with a momentum dependent fermion 
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Fig. 5. (a) Contribution to the 2PI effective action at two-loop order in the loop 
expansion or at next-to-leading order in the 1/Nf expansion, (b) The corresponding 
fermion self-energy and vacuum polarization at the same order. All internal lines 
denote the full propagators and all vertices are the bare ones. 

self-energy that violates the corresponding WT identity, is gauge independent 
up to corrections of higher order 

In Appendix B we briefly summarize some exact relations derived from the 2PI 
effective action that are useful to our discussion here. The interested reader 
is referred to the literature [63] for further details. For QED in a constant 
external magnetic field, the 2PI effective action is given by 

F[A,G,V] = S [A] - tTr log G- 1 - iTr S-\G - S) 

+ ^Tr logD- 1 + -TrD'^V - D) + T 2 [G,V], (5.3) 

where Sq[A] is the classical Maxwell action with A^ being the external gauge 
field, and the propagators are the same as those in Sec. 3. Specifically, we con- 
sider the LLLA hence the fermion propagators are those for the LLL fermions. 
In actual applications, the 2PI effective action is truncated at some order in 
a chosen expansion parameter. For the purpose of our discussion, it suffices 
to consider only the lowest nontrivial order truncation, namely, at two-loop 
order in the loop expansion in QED or at next-to-leading order [i.e., 0(iV2)] 
in the 1/Nf expansion in large- Nf QED. 

At the lowest nontrivial order in the loop expansion or in the 1 / Nf expansion 
T 2 [G, V] is given by 

V 2 [G,V] = ~J d A xd A ytT^G{x,y)YG{y,x)V^{x,y), (5.4) 
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where in the 1/Nf expansion e is understood as the rescaled coupling that 
remains finite and fixed in the limit Nf — > oo. Functional differentiations 
of T[A, G, T>\ with respect to the full propagators in the absence of sources 
leads to the SD equations for the respective propagators. The corresponding 
fermion self-energy and vacuum polarization are obtained by varying r 2 [G, T>\ 
with respect to the corresponding full propagators. The resultant SD equations 
are found to be exactly the same as those in the BVA with the corresponding 
fermion self-energy and vacuum polarization given by Eqs. (4.16) and (4.5), 
respectively. A diagrammatical representation of the 2PI effective action is 
depicted in Fig. 5, in which it is noted that the vertices that appear in r 2 [G, V\ 
as well as in the corresponding fermion self-energy and vacuum polarization 
are the bare ones. With the above observation one can conclude fairly easily 
that the 2PI effective action truncated at the lowest nontrivial order in the 
loop expansion or in the 1/Nf expansion resums identically the same infinite 
subset of diagrams as the SD equations truncated in the BVA. 

The direct correspondence with the 2PI effective action truncated at the low- 
est nontrivial order has the following important consequences. On the one 
hand, the WT identity in the BVA (4.4) guarantees that the truncation at 
the lowest nontrivial order in the loop expansion or in the 1/Nf expansion 
with a momentum independent fermion self-energy is a consistent truncation 
of the 2PI effective action in the momentum region relevant to chiral sym- 
metry breaking in a strong magnetic field. On the other hand, the facts that 
vertex corrections in the 2PI effective action start to appear at the next-to- 
lowest-nontrivial order in the loop expansion or in the 1/Nf expansion and 
that the truncated on-shell 2PI effective action (in general with a momentum 
dependent fermion self-energy) has a controlled gauge dependence with the 
explicit gauge dependent terms always appearing at higher order, provide an 
unequivocal justification of our argument that in a strong magnetic field the 
dynamical fermion mass to leading order in the gauge coupling or in the 1/Nf 
expansion can be reliably calculated in the BVA, and is gauge independent up 
to corrections of higher order. Whereas a direct verification by explicit calcu- 
lation is indispensable, to the best of our knowledge consistent truncations of 
the SD equations that include vertex corrections either within or beyond the 
LLLA have not appeared in the literature. Clearly such a task is beyond the 
scope of the current article, and will be the subject of further investigations. 



6 Summary and conclusions 

In this article we presented a critical and detailed study of chiral symmetry 
breaking in QED in a constant external magnetic field. The main goal is 
to determine in a gauge independent manner the dynamical fermion mass 
generated through chiral symmetry breaking. 
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We focused on the strong field limit as well as on the weak coupling regime. 
The former leads to a wide separation of energy scales such that the dynamics 
of fermion pairing is dominated by the lowest Landau level. The latter allows 
for a controlled, nonperturbative calculation of the dynamical fermion mass, in 
the sense that there exists a systematic expansion in powers of the gauge cou- 
pling (up to logarithms) such that contributions of leading order in the gauge 
coupling that arise form every order in the loop expansion can be consistently 
accounted for, while subleading contributions are suppressed by powers of the 
gauge coupling (up to logarithms). Consequently, in the weak coupling regime 
a gauge independent, consistent truncation of the SD equations is tantamount 
to a good and reliable approximation. 

The WT identity in the BVA is at the heart of our proof that the BVA is a 
consistent truncation of the SD equations within the LLLA. We first verified 
that, with the corrected orthonormal condition of the E p functions for the LLL 
fermions, the conclusion of Ref. [23] that in order to satisfy the WT identity 
in the BVA within the LLLA, the fermion self-energy has to be a momen- 
tum independent constant remains valid. Furthermore, we showed that the 
WT identity in the BVA can be satisfied only within the LLLA. The proof 
then proceeds with the assumption that the BVA is a consistent truncation 
such that the position of the fermion pole obtained therein is gauge indepen- 
dent. With this assumption, the WT identity in the BVA implies that the 
pole of the LLL fermion propagator is located at p| = — m 2 , with m being 
the momentum independent fermion self-energy as well as the physical, gauge 
independent dynamical fermion mass. The proof is completed by verifying 
that such an assumption does not lead to any inconsistency. This is achieved 
by showing that (i) the truncated vacuum polarization is transverse, (ii) the 
truncated fermion self-energy evaluated on the fermion mass shell is mani- 
festly gauge independent. In particular, we showed that the would-be gauge 
dependent contribution to the truncated fermion self-energy that arises from 
the gauge dependent term in the full photon propagator vanishes identically 
on the fermion mass shell. This detailed study leads to a gauge independent, 
nonperturbative calculation of the dynamical fermion mass to leading order 
in the gauge coupling. Furthermore, it allows one to identify unambiguously 
the infinite subset of diagrams that contribute to chiral symmetry breaking in 
a strong magnetic field. 

We made a detailed comparison between our results and those obtained in pre- 
vious works [18-20,22,24-29]. The gauge independent analysis we presented 
in this article shows clearly that the results as well as the conclusions of those 
earlier studies are inevitably gauge dependent and can be attributed to incon- 
sistent truncation schemes, gauge dependent artifacts, or both. Nevertheless, 
based on the gauge dependent results [24-28] obtained within the LLLA re- 
garding the reliability of a momentum independent fermion self-energy in the 
BVA in the momentum region \p^\ <C \eH\, we argued that in a strong mag- 
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netic field the dynamical fermion mass can be reliably calculated in the BVA 
and is gauge independent up to corrections of higher order. This is consistent 
with the underlying assumption of the LLLA that the kinematic region of low- 
energy photon exchange ?x <C \eH\) is precisely that which determines 
the dynamics of fermion pairing in a strong magnetic field. 

Motivated by the fact that the gauge dependence of the dynamical fermion 
mass calculated in the BVA but with a momentum dependent fermion self- 
energy is subleading, we established a direct contact between the SD equations 
truncated in the BVA and the 2PI effective action truncated at the lowest 
nontrivial order in the loop expansion as well as in the 1/Nf expansion. This 
direct correspondence, together with the facts that vertex corrections in the 
2PI effective action start to appear at the next-to-lowest-nontrivial order in 
the loop expansion or in the 1/Nf expansion and that the truncated on-shell 
2PI effective action has a controlled gauge dependence with the explicit gauge 
dependent terms always appearing at higher order, provides the justification 
of our argument that in a strong magnetic field the dynamical fermion mass 
can be reliably calculated in the BVA to leading order in the gauge coupling 
or in the 1/Nf expansion. 

In conclusion, the presence of a strong, constant external magnetic field pro- 
vides a consistent truncation of the Schwinger- Dyson equations in terms of 
the bare vertex approximation within the lowest Landau level approximation. 
This allows us to study in a gauge independent manner the physics of chi- 
ral symmetry breaking in QED in an external magnetic field and to obtain a 
solution for the physical dynamical fermion mass. The gauge independent ap- 
proach to the Schwinger-Dyson equations in gauge theories that we presented 
in this article is quite general in nature, and hence not specific to the prob- 
lem of chiral symmetry breaking in a strong magnetic field discussed here. 
We believe this approach will be useful in other areas of physics that also 
require a nonperturbative understanding of gauge theories, such as the stud- 
ies of hadronic structure and the connection between quark confinement and 
chiral symmetry breaking in nuclear and particle physics [50,64], the physics 
of superconductivity and superfluidity in condensed matter systems [65,66], 
as well as the equilibrium and nonequilibrium properties of hot and dense 
matter [67,68]. 
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A Vacuum current in a constant magnetic field 



Using the Ritus E p functions and the general structure of the fermion self- 
energy, we shall show in a nonperturbative manner that the vacuum current 
vanishes identically in a constant external magnetic field. The vacuum current 
is defined to be 

J» = ietr^G(x,x) 

where in the second equality we have used the momentum representation of the 
full fermion propagator given by Eq. (2.32). It is noted that in this appendix 
the fermion propagator and the fermion self-energy will not be restricted to 
those for the LLL fermions. 

With the properties (2.14) and (2.15), the integration over p 2 , one of the spin 
sums, as well as the trace can be done easily, leading to 



d 2 Pll [i + An(p)]pf 



(A.2) 

where p?, = —pi +p\,p\ = 2\eH\l, and use has been made of the Dirac matrix 
structure of the fermion self-energy in Eq. (3.13). Since A\\ t ±(p) and B{p) are 
functions of the longitudinal and transverse momentum squared, one finds 
that vanishes identically. The same conclusion also applies to the vacuum 
current in the LLLA, in which the sum over I is restricted to the / = term 
and that over a is restricted to a — sgn(cff). 

We note that the vanishing of the vacuum current in a constant external 
magnetic field is a consequence of the Lorentz boost invariance along, as well as 
the rotational invariance about, the direction of the constant external magnetic 
field, namely, the ^-direction. 
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B 2PI effective action 



In this appendix we briefly summarize some exact relations derived from the 
2PI (two-particle-irreducible) effective action [63] that are useful to our dis- 
cussion. For presentational simplicity we consider here only bosonic fields, the 
generalization to fermionic fields is straightforward. The generating functional 
for connected correlation functions is defined as 

Z[J,K] = e iW ^ 

= J T> (pe i ( s l<p}+ J i<P i H<P iK ii<P j ) : (B.l) 

where S[tp] is the classical action, tp l denotes bosonic fields and J (K) is 
the auxiliary local (bilocal) source. We use a shorthand notation in which 
Latin indices stand for space-time variables as well as internal indices, and 
summation and integration over repeated indices are understood, e.g., 

JW* = J d A xJ(x)<p(x). (B.2) 



The mean field 0* = (cp 1 ) and the connected two-point function G' lj = (T(p' l (p j ) — 
(ip l )(ip j ) are obtained by functional differentiations of W[J, K] with respect to 
the local source J as 

0{lJi) d(lJi)d(lJj) 

Functional differentiations of W[J, K\ with respect to the bilocal source K may 
generate also disconnected correlation functions. For example, differentiating 
once with respect to K leads to 

5W 1 

'sm7) = 2 m ' +G " y (R4) 

A functional Legendre transform in the mean field <fi and the connected two- 
point function G leads to the 2PI effective action 

G] = W[J, K] - U l - \k % 0<\? + G ij ). (B.5) 

From Eqs. (B.3) and (B.4) one can derive the relations 

6(f) 1 3 5G % 3 2 3 y J 

The 2PI effective action can be cast into a very convenient form which has a 
simple diagrammatical interpretation in terms of 2PI diagrams (for fermionic 
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fields the factors of 1/2 before the trace terms are replaced by —1) 

T[<j>, G] = S [<f>] + l -Tr logG" 1 + '-Tr G \G - G ) + r 2 [0, G], (B.7) 

where So is the free part of the classical action and Go is the bare two-point 
function 

The functional r 2 [0, G\ is the sum of all 2PI skeleton vacuum diagrams with 
bare vertices and full propagators. Here, skeleton diagrams are those without 
self-energy insertions. 

The equations of motion for the mean field <fi and the connected two-point 
function G are determined by the stationary condition, i.e., from the implicit 
functional equation (B.6) for vanishing sources J and K, as 

^M = , ^Mi=0. (B.9) 

5(p l SG % 3 v ' 

With Eq. (B.7) the equation of motion for the connected two-point function 
in Eq. (B.9) becomes 

G'j 1 = G^) + Eij, (B.10) 

where £ is the one-particle-irreducible self-energy (for fermionic fields the 
factor of —2 is replaced by 1) 

pun 

It is noted that Eq. (B.10) is the Schwinger-Dyson equation for the full prop- 
agator G. 
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